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ABSTRACT 

We present the analysis of 4.5 years of nearly continuous observations of the classical Cepheid 
Polaris, which comprise the most precise data available for this star. We have made spectroscopic 
measurements from ground and photometric measurements from the WIRE star tracker and the SMEI 
instrument on the CORIOLIS satellite. Measurements of the amplitude of the dominant oscillation 
(P = 4d), that go back more than a century, show a decrease from Ay = 120 mmag to 30 mmag 
around the turn of the millennium. It has been speculated that the reason for the decrease in amplitude 
is the evolution of Polaris towards the edge of the instability strip. However, our new data reveal an 
increase in the amplitude by ~ 30% from 2003-2006. It now appears that the amplitude change is 
cyclic rather than monotonic, and most likely the result of a pulsation phenomenon. In addition, 
previous radial velocity campaigns have claimed the detection of long-period variation in Polaris 
(P > 40 d). Our radial velocity data are more precise than previous datasets, and we find no evidence 
for additional variation for periods in the range 3-50 d with an upper limit of 100ms . However, in 
the wire data we find evidence of variation on time-scales of 2-6 days, which we interpret as being 
due to granulation. 

Subject headings: stars: individual: a UMi (HD 8890; Polaris) - stars: variables: Cepheids 



1. INTRODUCTION 

In addition to being arguably the most famous and 
in practice useful star other than the Sun, Polaris has 
a number of properties that may provide insights that 
are important to stellar astrophysics. It is the near- 
est and brightest classical Cepheid, oscillating in a sin- 
gle mode of pulsation with a period around four days. 
The H IPPARCOS parallax constra ins its luminosity and al- 
lowed [IiasF&ICa|chi)oIl ()1997j ) to argue that the mode 
of pulsation must be the first overtone, which is up- 
held by the recent reeval uation of the HIP PARCOS data 
van Leeuwen et aLl 120071 ). As discussed by lEvans et al.l 



TABLE 1 

Observing log. The point-to-point scatter, a, is given in 

MMAG (WIRE + SMEI PHOTOMETRY) AND KMS -1 (SPECTROSCOPY). 











Data 




Source 


Date start 


Date end 


T obs [d] 


points 




WIRE 


2004 Jan 17 


2004 Feb 15 


29.0 


41 162 


0.88 


WIRE 


2004 Jul 10 


2004 Jul 30 


20.2 


33 937 


0.30 


WIRE 


2005 Jan 31 


2005 Feb 12 


12.6 


14 240 


0.29 


SMEI 


2003 Apr 6 


2006 Dec 31 


1365.4 


13 543 


6.4 


AST 


2003 Dec 25 


2007 Oct 26 


1401.0 


517 


0.15 


D89 


1987 May 6 


1988 Jan 3 


241.0 


175 


0.44 


H00 


1991 Nov 21 


1993 Aug 2 


640.4 


42 


0.07 



2002) , Polaris has a number of unusual pulsation proper- 



ties, including a very small pulsation amplitude, and, as 
for other overtone pulsators, it has a rapid peri od change. 
It has bee n clear from stud ies as far back as iParenagol 
(|1956ft and ISzabadosl (|1977f ) that overtone pulsators ( "s 
Cepheids" as they were then called) had rapid period 
changes, more rapid than can be explained by evolution 
during a second or third crossing of the i nstability strip. 

It was found by lArellano Ferrol (|1983l ) that the main 
period of Polaris is increasing (316 s per century) and 
that its peak-to-peak amplitude has decreased signifi- 
cantly, from about 140 to 70 mmag (Johnson B filter), 
based on photoe lectric measurements col lected in the pe- 
riod 1930-1980. iKamper k, Fernid (|1998l) analysed radial 
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velocity measurements from 1900 to 1998 and could con- 
firm the decrease in amplitude. However, their radial 
velocity data from the end of the pe riod showed t h at the 
decrease had apparently stopped. iTurner et al.l (2005) 
analysed both radial velocity and photometry data and 
found evidence for a sudden change in oscillation period 
around 1963-1966, and that the amplitude change be- 
came steeper at the same time. 

Our aim is to shed light on the properties of the oscil- 
lation of Polaris, based on the analysis of a nearly con- 
tinuous dataset spanning 4.5 years. We have photome- 
try from two satellite missions and simultaneous spectro- 
scopic monitoring from the ground, each dataset being 
superior to previous data available for Polaris. 

2. OBSERVATIONS AND DATA REDUCTION 

Polaris (a Ursae Minoris) was observed with the star 
tracker on t he Wide-field InfraRed Explorer (wire; 
iBrunttJ I2007T) satellite in January and February 2004, 
July 2004, and February 2005. These runs lasted about 
4, 3 and 2 weeks, respectively. In addition, Polaris was 
monitored using 3.8 years of nearly continuous photom- 
etry from the SMEI instrument on the CORIOLIS satel- 
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Fig. 1. — The photometry datasets from SMEI and WIRE and the radial velocities from AST. The bottom panel is an expanded view of 
40 days of observations surrounding the last run with WIRE. The thick gray curve is the fit of a single sinusoid to the SMEI data. 



TABLE 2 

Times of maximum light in the smei and wire datasets. 



t E 


Epoch 


Source 


[HJD -2, 453, 000] 






391.327 


6329 


SMEI 


395.339 


6330 


SMEI 


399.387 


6331 


SMEI 


403.323 


6332 


WIRE 


407.274 


6333 


WIRE 


411.262 


6334 


WIRE 


415.418 


6335 


SMEI 


419.160 


6336 


SMEI 


423.215 


6337 


SMEI 


427.152 


6338 


SMEI 



lite. These observations were obtained between April 
2003 and the end of 2006. We further used the 2 m Ten- 
nessee State University Automatic Spectrosc opic Tele- 
scope (AST; lEaton fc Williamsonll200l l200l to collect 
517 high dispersion spectra of Polaris over a period of 
3.8 years, from late 2003 to late 2007. A log of the pho- 
tometric and radial velocity observations is given in Ta- 
ble [T] We also list datasets from two previous radial 
velocity campaigns that we ha ve used for compariso n in 
our analysis in Sect. PS (D89: iDinshaw et al.lTl989l and 
H00: lHatzes fc Coch7anlf2000l) . 

The complete photometric light curve from SMEI and 
WIRE and the radial velocity data from AST are shown 
in the top panel in Fig. [T] Note that different units for 
the photometry (mmag) and velocities (kms" 1 ) are given 
on the left and right abscissa, respectively. The right ab- 
scissa is adjusted by the ratio of the measured amplitudes 
in the AST spectroscopy and the smei photometry. The 
bottom panel in Fig. [1] shows the details of the variation 
during 40 days. It shows the last run done with WIRE 
and illustrates the typical coverage with SMEI and the 
AST. The thick gray curve is the fit of a single sinusoid 
to the complete SMEI dataset. 

The wire dataset consists of about three million CCD 
stamps extracted from the 512x512 CCD SITe star 



tracker camera. Each window is 8x8 pixels and we 
carried ou t aperture ph o tomet ry using the pipeline de- 
scribed bv lBruntt et al.l (|2005h . The resulting point-to- 
point scatter range from 0.3 to 0.9 mmag (see Tabled]). 
The noise is somewhat higher than the Poisson noise, e.g. 
the high noise in the first WIRE dataset is due to a high 
background sky level. 

Data from the AST consist of echelle spectra cover- 
ing the wavelength range 5 000-7100 A at a resolu- 
tion of about 30 000. The velocities are derived from 
the correlation between the observed spectrum and a 
list of 74 mostly Fei lines. More details about the 
pip eline used to extract the r adial velocities are given 
by lEaton fc Williamson! (|2007l ). We find that the ve- 
locities deviate slightly from the IAU velocity system 
(At) = -0.35 ± 0.09km s" 1 ). The drift of the velocities 
due t o the motion of Polaris in its binary orbit (|Kamperl 
Gill was subtracted before the time-series analysis was 
carried out. 

Th e Solar Mass E j ection Imager (smei) (|Evles et al.l 
2003^ Uackson et all 12001 is a set of three cameras 
mounted on the side of the CORIOLIS satellite. Each cam- 
era covers a 3 x 60 degree strip of sky and this is projected 
onto 1260 x 40 pixels of a CCD. As the nadir-pointing 
satellite orbits every 101 minutes in its Sun-Synchronous 
polar orbit the cameras take continuous 4-second expo- 
sures, which result in roughly 4,500 images per orbit cov- 
ering most of the sky. After bias and dark removal and 
flat-fielding processing, the pixels about the target star 
are selected from about ten images that contain the star 
for each orbit. These pixels are then aligned and a stan- 
dard PSF is fitted by least-squares. Thus, the result is 
one brightness measure every 101 minutes. There are 
problems with particle hits and a complex PSF, which 
shows temporal variations that are not fully understood. 
The short-term (t < 1.0 d) accuracies are of the order of 
a few mmag, but there are longer term systematic errors 
of the order of 10 mmag on the timescales of days and 
months, which are not yet well understood. Allowing for 
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Year: 1940 1960 19B0 2000 




500 1000 




-30000 -20000 -10000 

Time = HJD - 2,453,000 [days] 

Fig. 2. — O — C diagram. The solid line is a fit to the data taken 
after 1965. The inset shows the details of the new O — C data from 
SMEI and WIRE. 

time allocated for calibrations and satellite operation, a 
time coverage of about 85% was maintained. 

3. PERIOD CHANGE OF THE MAIN MODE 

Observed epochs of maximum light in Polaris, span- 
ning more than a century, have indicated a significant 
change in period (I Arellano Ferrolll983l : lFernie et alJT993t 
iTurner et al.ll2005[ ). These studies have considered O — C 
diagrams 7 and fitted a parabola, which is equivalent to 
assuming the period changes linearly with time. We 
have combin e d new measured epochs with those from 
ITurner et al.l |2005). A sudden change in the period has 
been noted in the 1960's and therefore we have only used 
epochs observed since 1965. 

We measured 223 and 16 epochs of maximum light 
in the SMEI and WI RE datasets, respectiv ely. From the 
O — C analysis by lArellano Ferrol ()1983f ) we expected 
the period to decrease by 14 seconds over the 4.5 years 
time span of our dataset. Thus, when predicting the 
epochs of maximum light we could safely assume a con- 
stant period. The epochs were predicted by fitting a 
single sinusoid to the SMEI data. The fit is of the form 
S(t) = Asin[2?r(/ (t — to) + (f>)], where A is the ampli- 
tude 8 , t is the time (Heliocentric Julian Date), <o is the 
zero point, P = l/f is the period, and (j> is the phase. 
The result is P = 3.972111 ± 0.000054 days and phase 
(/) = 0.2158 ± 0.0022 for the zero point t = 2,530,000. 
This fit is shown with a thick gray line in the bottom 
panel in Fig. [1] The uncertainties were determined from 
simulations as discussed in Sect. 14.11 

At each epoch, E, we selected all data points within 
half a period (\t — E * P\ < P/2). For the SMEI data we 
required that at least 30 data points were available with 
at least 10 data points before and after the maximum. 
To estimate the time of maximum light, we fitted these 
data with a sinusoid with fixed frequency and amplitude, 
while the phase was a free parameter. In Table [2] we list 
these times, including the epoch number following the 
definition by ITurner et al.l (|2005l ). and the source of the 
data. The complete table is available in the on-line ver- 
sion of the paper, while the times listed here correspond 
to the time interval covered in the bottom panel in Fig. [TJ 

Since the quality of the three datasets are very differ- 
ent, we computed weights based on the uncertainty of 

7 O — C diagrams display observed times of maximum light minus 
calculated times, assuming a constant period, plotted vs. time. 

8 Previous studies on Polaris have used peak-to-peak amplitudes. 
We use amplitudes unless otherwise specified. 



TABLE 3 

Results of the frequency analysis for the five datasets. 



Source 


/ May- 1 ] 


a [mmag/kms ] 


[0..1] 


WIRE 


0.251732(14) 


12.26(18) 


0.2195(34) 


SMEI 


0.2517553(34) 


12.63(13) 


0.2158(22) 


AST 


0.2517506(43) 


0.998(14) 


0.1550(34) 


D89 


0.25174(14) 


0.808(58) 




H00 


0.251722(22) 


0.781(22) 





the epoch times, tg- This was determined by calculating 
the point-to-point scatter after subtracting a parabola 
fitted to the O — C data of each dat aset. The uncer- 
tainty on t E from ITurner et alj (|2005l ) is a(t E ) = 0.18d 
while for SMEI and wire data the uncertainties are 0.12 d 
and 0.014 d. Finally, we calculate the O — C values 
using weights, l/o~ 2 (tE), and the most accurate value 
of the period, which is the one from SMEI: O — C = 
(tE — to) — (E — E ) ■ -Psmei- We chose the refer- 
ence epoch to be Eq = 6333, and using the fit to the 
SMEI data found above, this corresponds to the time 
t = 2, 453, 407.277 ±0.009. 

Our final O — C diagram is shown in Fig. [2j The recent 
data continue to match the smooth period change as in- 
dicated by the weighted parabolic fit, corresponding to 
a period increase of dP/dt = 365 ± 27 sec/century. We 
note that the most recent SMEI data lie systematically be- 
low the fit, which is tightly constrained by the extremely 
well-defined times from WIRE. Note that if we only use 
data more recent than 1985, we conclude that the period 
has not changed, i.e. dP/dt = — 19 ± 95 sec/century. 

4. AMPLITUDE CHANGE OF THE MAIN MODE 

To analyse the time series in the frequency domain, 
we calculated the Fourier amplitude spectra for each 
dataset. In the top panels in Fig. [3] we show the ampli- 
tude spectra for the photometry from the wire (red) and 
SMEI (green) and the radial velocities from AST (black). 
The inset shows the prominent peak around 0.25 c day -1 , 
corresponding to the known 4-day period. Due to the 
long gaps between the three wire datasets the spectral 
window shows a more complex pattern than the nearly 
continuous datasets from SMEI and AST. 

We fitted a single sinusoid to the dominant peak and 
the results for each dataset are given in Tabled We list 
the phase of the WIRE, SMEI and AST datasets, relative to 
the zero point to = 2, 453, 000. It is seen that the frequen- 
cies are in very good agreement. Although the point-to- 
point precision of the wire data is superior, the long 
span of the SMEI and AST datasets means the frequency 
is determined more accurately by a factor of 3-4. We also 
note that the phases from the fit to the SMEI and wire 
photometry are in good agreement, but the phase of the 
fit to the AST data is quite different. This corresponds 
to a shift in the times of maximum in the flux and radial 
velocity: HJD( F maY ) = HJD(.RV maY ) + (1 .744 ± 0.016 d). 
In comparison, Mos kalik fe Qglozai ([2000) found the off- 
set to be 1.682 ± 0.017 d u sing combined pho t omet ry and 
radial velocity data from lKamper fe Fern ie (1998). The 
shifts appear to be mar ginally different (2 . 7 a), b ut we 
note that the combined iKamper fc Ferniel (|1998f ) data 
is calibrated to Johnson V, while SMEI has a filter re- 
sponse roughly cent ered at Johnson R but being wider 
(|Tarrant et alj|2007D . 

We have subtracted the main oscillation mode and the 
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Fig. 3. — Amplitude spectra of the AST radial velocity data and the photometric data from SMEI and WIRE. The top panel is for the raw 
data and the bottom panel is after subtracting the main mode at ~ 0.25 c day -1 . The insets show details of the amplitude spectra. 

. Year: . ao , M !°°! z ™ 2 ^Z 30 pulsation cycles, and the three wire datasets were 

Iia analysed independently. Each subset is fitted by a sin- 

- > gle sinusoid and we analysed the output parameters. We 

__,-A-x.s ^ made the analysis for two different assumptions: 




500 1000 
HJD - 2,453.000 [days] 



Fig. 4. — Increase in the amplitude of the 4d main mode mea- 
sured over 4.5 years with WIRE (small dots), SMEI (x), and AST (•). 



residual amplitude spectra are shown in the bottom panel 
in Fig. [3l All spectra show a higher amplitude towards 
lower frequencies. This could be either due to instru- 
mental drift or long-period variation intrinsic to Polaris, 
which we discuss in Sect. [SJ The left inset shows that 
the peaks in the SMEI data do not coincide with peaks 
in the AST data, so we cannot claim they are due to 
coherent pulsations. The right inset shows the details 
around 0.25 c day -1 , where two significant residual peaks 
are seen in the SMEI and AST datasets. They are an in- 
dication that the frequency, amplitude or phase of the 
oscillation has changed during the ~ 4 years of observa- 
tion. 

To test the possible change in either frequency, ampli- 
tude or phase of the main mode, we divided the datasets 
into subsets. Each of the SMEI and AST subsets contain 



• We assumed that the frequency, amplitude and 
phase change with time, so each is a free parameter. 

• We assumed that only the amplitude changes with 
time, so the frequency and phase were held fixed. 

The first assumption allowed us to test the stability 
of the mode. We found that within the uncertainties 
both the frequency and phase do not change over the 
time span of the observations. However, the amplitude 
increases monotonically, and this result is confirmed by 
the analysis under the second assumption. A comparison 
of the measured amplitudes under the two assumptions 
give nearly identical amplitudes for all subsets. 

In Fig. [4] we show the result obtained under the second 
assumption, i.e. for fixed frequency and phase. The x 
symbols are the SMEI data, the • symbols are the AST 
data, and small dots are for the three WIRE runs. The 
dashed line is a weighted linear fit to the SMEI data 
yielding the amplitude change in mmag: j4 p hot(t) = 
(10.94 ± 0.08) + (2.56 ± 0.13) • 10" 3 (t - to), where t is 
the HJD with zero point to = 2,453,000. The weighted 
fit to the AST data (solid line) gives: A RV (t) — (0.90 ± 
0.01) + (1.45 ± 0. 15) ■ 10~ 4 (t - 1 ) km s" 1 . From this anal- 
ysis we find that from 2003-7 the amplitude measured 
in flux and radial velocity has increased by 34 ± 2% and 
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Fig. 5. — Peak-to-peak radial velocit y amplitude for the mai n 
4d mode. The solid line is the fit from Kampcr & Fcrnie (1998). 
The inset shows the most recent dat a in detail. Our mea sure- 
ments of the ampl i tude based on the Dinshaw ct al. (198||) and 
IHatzes fc Cochran (20Q3) are marked by arrows. 

24 ± 3%, respectively. The single mode in Polaris has 
very low amplitude and from linear theory it is expected 
that the rate of increase is the same in photometry and 
radial velocity. The measured increase over four years is 
only marginally different (8.1 ± 3.1% or 2.6 ct), and con- 
tinued monitoring is required to confirm this tentative 
result. 

In Fig. O the measured peak-to-peak amplitudes in ra- 
dial velocit y from AST are compared with results from the 
literature. iKamper fc Ferniel (|1998h found a monotonic 
decrease in the amplitude over the past 100 years 9 To 
avoid confusion about the conversion factor, we compare 
our results directly in radial velocity, and we have repro- 
duced their fit as the solid black line. The dashed line 
is an extrapolation, and as noted by IKamper fc Ferniel 
(1998), zero amplitude would be predicted in the year 
2007. However, from four yea rs of monitoring (1 9 93.15- 
1996.96) in radial velocity, IKamper fc Ferniel (|1998[ ) 
found that the peak-to-peak amplitude was nearly con- 
stant. This is in agreement with the amplitudes we 
fin d from the an a lysis (see Sect. 15.11) of the dataset s 
bv lDinshaw et alj (|1989D and IHatzes fc Cochran! (|2000f ). 
From the most recent high-precision data, it is evident 
that the amplitude in Polaris was constant in the period 
1987-1997, while the increase we report here marks a 
new era in the evolution of Polaris. 

4.1. Uncertainties in frequency, amplitude and phase 

The uncertainties in the frequency, amplitude and 
phase in Table [3] are determined from realistic simu- 
lations of each dataset. The times of observation are 
taken from the observations and the simulations take 
into account the increase in noise towards low frequen- 
cies. We used the approach described by iBruntt et al.1 
(2007) (their Appendix B) to make 100 simulations of 
each dataset, and the uncertainties are the rms value on 
the extracted frequency, amplitude and phase. 

We made simulations of the wire, SMEl and AST 
datasets and also the radial velocity datasets from 

9 IKam per & Fcrnie ( 1998) converted their radial velocity am- 
plitudes to photometric amplitudes using an empirical factor 
50 km s^ 1 mag -1 . The ratio of the amplitudes measured in our 
radial velocity and the SMEl photometric data is Aptv/^phot = 
79. 0± 1.4kms _1 . The reason for the value being high er is that the 
filter r esponse of SMEl is roughly Johnson R while Kampcr fc Ferniel 
(1998) used Johnson V. 




0.00 0.10 0.20 0.30 

Frequency [c/day] 

Fig. 6. — Residual amplitude spectra of the three radial velocity 
datasets from AST (black), Hatzes & Cochran (2000; blue), and 
Dinshaw et al. (1989; red). The low frequency peaks claimed in 
the previous studies are marked by arrows. 

iDinshaw etaTl (|1989l ) and IHatzes fc Cochran! (|2000l) . 
which we use for comparison in Sect. 15.11 In all 
cases the uncertainties are somewhat higher than the- 
oretical estimates that assume th e noise to be white 
([Montgomery fc O'Donoghuel 119991 ). This is especially 
the case for the wire photometry due to low-amplitude 
variation on timescales comparable to the 4 d period, as 
will be discussed in Sect. 15.21 

5. EVIDENCE FOR ADDITIONAL INTRINSIC VARIATION 

In the following we will look at the evidence for varia- 
tion in Polaris beyond the 4d main mode. The analysis 
is based on the AST and wire datasets and data from two 
previous radial velocity campaigns. The AST dataset is 
the most accurate at low frequencies, and can be used 
to look for long periods (Sect. 15. 1|) . The wire dataset 
consists of three week-long runs with very high preci- 
sion and can be used to study low-amplitude variation 
(Sect. 15. 2[) . The noise in the SMEl dataset is distinctly 
non- white, as mentioned briefly in Sect.[2j Therefore, we 
did not use this dataset in the following investigations, 
except for the subtraction of the 4 d mode, where SMEl 
provides the most accurate frequency and phase. 

5.1. Long-period variation 

In addition to the 4 d mode, a few studie s have claimed 
varia tion at long pe riods, e.g. P = 9 . 75 d (IKamper et al.1 
fl98l . 45.3 ± 0.2 d IDinshaw et all (fl989T) and 40.2 ± 
0.7 d and 17.2 d (|Hatzes fc Cochran! 12000) . To con- 
fir m these claims we r eanal ysed the original dataset s 
of IDinshaw et all (|1989f l and IHatzes fc Cochran! <|2000f ). 
The basic properties of these datasets are listed in Ta- 
ble [Q 

In Fig. [6] we compare the amplitude spe c tra o f the 
spectroscopic data from A ST, IDinshaw et al.1 (|1989f l and 
IHatzes fc Cochranl (|2000f ). after subtracting the 4d pe- 
riod and the long-period trend due to the binary orbit. 
The AST spectrum was calculated taking the amplitude 
increase into account (cf. Sect. 01 . As a result, the resid- 
ual double peak seen at ~ 0.252 cday -1 is no longer vis- 
ible (compare Fig. [6] with the bottom panel in Fig. [3]). In 
the AST amplitude spectrum we see no significant peaks 
from 0.02-0.3 cday -1 (P = 3-50 days). We set an up- 
per limit on the amplitude at 100ms -1 , which is four 
times the average level in the amplitude spectrum in this 
frequency interval. 

The dataset by IHatzes fc Cochranl (|2000D consists of 
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Time = HJD - 2,453,000 [days] 

Fig. 7. — Residuals in the three WIRE light curves of Polaris 
and the comparison star HD 166205 (large • symbols offset by 
3mmag). The dashed curve is a sinusoid with the 4d period and 
an amplitude set arbitrarily to 2 mmag. 

42 data points distributed unevenly over 640 days, pro- 
viding a very complicated spectral window. The dataset 
comprises very precise radial velocities collected with an 
iodine cell as a re f erence (rms residuals are 70 ms -1 ). 
lHatzes fc Cochran! (|2000f ) detected two peaks with al- 
most equal amplitude at /^ 00 = 0.02 49 and f^ 00 = 
0.0581 cday -1 . In the amplitude spectrum in Fig. [6] we 
find only one of these peaks to be significant (/goo = 
0.05834 ± 0.00014 cday" 1 with amplitude A = 163 ± 
25ms _1 ). We made 100 simulations of the data includ- 
ing this frequency, white noise, and a 1// noise compo- 
nent consistent with the observations. The inserted /h 00 
frequency was only recovered in 49 of the 100 simula- 
tions. Interestingly, in the other simulations the highest 
peaks were clustered close to /^ 00 (0.025 ± 0.007c day -1 
with amplitude A = 1 67±22ms~ 1 ), indi c ating that the 
two peaks detected by lHatzes fc Cochran! (|2000[) are due 
to the complicated spectral window. 

The dataset by IDinshaw etall ([1989) comprises 175 
velocities distributed over 241 days. The precision is 
440 ms -1 , as estimated from the rms of the residuals 
after subtracting the 4d mode. In addition to this mode 
we detect a peak at 1. 02199±0. 00024c day" 1 with strong 
aliases at ±1 c day -1 . IDinshaw et al.l ( 1989) found one of 
the alias peaks (/ D8 9 = 0.0221±0.0001cday" 1 ) to be the 
highest and suggested that it was intrinsic to Polaris. We 
measure the amplitude to be A = 0.542 ± 0.067 km s _1 , 
but such a strong signal is not seen in the more recent 
dataset by lHatzes fc Cochran! (j2000l ) or in any of our 
datasets (see Figs. [3] and EJ). For these reasons we be- 
lieve that this peak is a leday -1 artifact, likely caused 
by instrumental drift in combination with the spectral 
window. 

5.2. Low-amplitude variation 

The residuals in the wire light curves, after having 
subtracted the 4 d period, are shown in Fig. [7l The drifts 
seen in the wire light curves correspond to periods from 
2-6 d, i.e. similar to the 4d main mode. Before interpret- 
ing this variation, it is important to investigate if these 
variations are due to improper subtraction of the main 
period or due to instrumental drift. 

To test the first caveat, we compared the residuals 
when subtracting the fit to the SMEl data and when sub- 
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Fig. 8. — Power density spectrum (PDS) of the residual light 
curves from WIRE shown in Fig. [jj The dashed line is a scaling 
of the granulation measured in the Sun. The arrows mark the 4 d 
main period and frequencies that are due to the WIRE orbit. 

tracting fits to the individual wire datasets. We find 
that the residuals are quite similar, and the conclusions 
reached here are not affected by the adopted approach. 
In the following, we have used the most accurate value 
for the frequency and phase, which is from the fit of the 
main mode to the SMEl data, while we used the ampli- 
tudes fitted to the individual wire datasets. 

The second caveat is instrumental drift and we tested 
this by using a comparison star. During the observations 
with the WIRE star tracker four other stars were moni- 
tored on the same CCD. However, not all are suited as 
comparison stars: HD 5848 is a bright (V = 4.3) K giant 
star clearly showing solar-like oscillations f Stel lo et alJ 
[20081 ). HD 51802 and HD 174878 are faint (V = 5.1 and 
6.6) M giants showing variation with long periods. The 
only suitable comparison star is the Al V star HD 166205 
(S UMi; V = 4.4). In Fig.[7]we also show the light curve 
of this star with large • symbols. The star is two mag- 
nitudes fainter than Polaris, so we binned the data col- 
lected within each orbit (P or bit — 93min) to be able to 
see any low-amplitude variation. 

The clear variation in the residuals of the wire 
photometry is not seen in the hot comparison star 
HD 166205. The fact that it is seen in three light curves 
obtained with wire at three different epochs spanning 
one year gives us some confidence that the signal is in- 
trinsic to Polaris. We will discuss two possible explana- 
tions for the observed variation here, namely granulation 
and star spots. 

To investigate in detail how the variation in the wire 
light curves depends on frequency, we show their power 
density spectra (PDS) in Fig. [5] in a logarithmic plot. 
The virtue of the PDS is that one can directly compare 
the properties of datasets, which have different temporal 
coverage and noise characteristics. Each PDS shows a 
clear increase from the white noise level around lOmHz 
towards low frequencies. The arrows mark the harmonics 
of the orbital frequency of wire (/wire = 178 /iHz), one 
third of the WIRE orbital frequency, and 2 c/day. These 
frequencies are observed in almost all wire datasets, 
and are due to a combination of a low duty cycle (typi- 
cally 20-40%) and scattered light from earth shine. Also 
marked is the location of the 4d mode (/ = 2.91 /iHz), 
which has been subtracted. 

To see if granulation could explain the increase in 
the PDS towards low frequencies, we have used a scal- 
ing of the granulation observed in the Sun. This is 
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based on VIRGO satellite observations of the Sun-as-a- 
star, and the scaling is done both in terms of amplitude 
and timescal e, following Kjelds en & Bedding (in prepara- 
tion; see also lStello et al.f2007f h This prediction is shown 
as the dashed line in Fig. [8] This is a scaling over three 
orders of magnitude in luminosity from the Sun to the 
super giant Polaris, so it is intriguing that the prediction 
of the granulation signal agrees with the observations 

within a factor of 3.0 ± .5 in amplit ude. 

Bot h iDinshaw et al.l (|1989f ) and lHatzes fc Cochran! 
(2000) discussed evidence for star spots on Polaris. We 
argue that this is not the cause of the observed variation 
in the wire data, since the time-scale is only a few days. 
That would imply very rapid rotation, which is unlikely 
for such an evolved star. 

6. DISCUSSION AND CONCLUSION 

From 4.5 years of intensive monitoring of Polaris we 
find that the amplitude of the 4-day main mode has in- 
creased steadily by about 30% in both radial velocity 
and flux amplitude. The rate of increase in the amplitude 
from 2003 -2007 is slightly st e eper t han the decrease from 
the fit bv lKamper fe Fernid (|1998( 1 for the period 1980- 
1994. Other sources hav e also found a rece nt increase in 
the amplitude of Polaris (|Engle et al.ll2004t Turner 2007, 
private communication) . 

The result that the amplitude of Polaris is now in- 
creasing has implications for earlier explanations of the 
change in amplitude as a cessation of pulsation due to 
the Cepheid's evolution towards the edge of the instabil- 
ity strip. At the very least, whatever the process is, it 
is not a simple monotonic progression through the HR 
diagram, unless this stage of evolution is more complex 
than previously thought. The recovery of the amplitude 
suggests that the phenomenon is cyclic. As such, it is 
likely to be associated in some way with pulsation rather 
than with evolution. 

A possible explanation for the increase in ampli- 
tude could be the beating of two closely spaced modes 
(jBreger fc Kolenbergl [2006) . To exhibit a beat period 
as long as the amplitude variation of Polaris (likely a 
few 100 years), the modes would have to be very closely 
spaced indeed. Among classical Cepheids, amplitude 
variation is extremely uncommon. The onl y case where 
it is firmly established is in V 473 Lyrae (jBurki et al.l 
1986), where the variation occurs over only a few years. 
In RR Lyrae stars, the Blazhko effect is well known, and 
thought to be the result of mode beating, although it is 
not completely understood. Since the amplitude changes 
in Polaris are well established, they require further ob- 
servations and consideration theoretically to unravel the 
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cause. 

There are other interesting aspects to the main oscil- 
lation mode. The characterization of the period change 
is puzzling since analysis of O — C diagrams spanning 
more than a century reveal that the period change is 
not linear. There may have been a "glitch" in both 
the period change a nd the amplitude in the mid-1960's 
(|Turner et a 1.1 [2005 h rather than smooth changes. We 
measured 239 epochs of maximum light but the time 
span of the observations of 4.5 years is too short to inves- 
tigate the period change. This is another aspect of the 
pulsation that demands further observation. 

A few radial velocity campaigns have claimed the pres- 
ence of additional long-period variation in Polaris. We 
have analy sed the original da t asets by Dinsha w et al.l 
(1989) and lHatzes fc Cochranl (|2000f ) and we conclude 
that these long periods are li kely spurious detecti ons 
caused by instrumental drifts (|Dinshaw et al.l IT989I ) or 
insuf ficient data leading to a complicated spectral win- 
dow (jHatzes fc Cochranl[200Q] ). From our 3.8 yr of mon- 
itoring with AST, we set an upper limit on the variation 
in radial velocity at 100 ms -1 for periods in the range 
3-50 d (except for the 4 d main mode) . 

In the WIRE data we find evidence of low-amplitude 
variation (peak-to-peak 2mmag) at time scales of 2-6 
days, which are likely intrinsic to Polaris. We applied 
a simple scaling of observed solar values for the charac- 
teristic timescale and amplitude of the granulation. Al- 
though this is a scaling over three orders of magnitude 
in luminosity, the prediction agrees with the observed 
variation in Polaris within a factor 3.0 ± 0.5. Thus, we 
conclude that the variation in the wire data could be 
due to granulation. 
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